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I Introduction 



There is a long-held belief that in quantum theory of gravitation space-time must 
change its nature at distances comparable to the Planck scale. In order to model 
such a situation one can invoke the Heisenberg uncertainty rules. In their standard 
form they make the notion of classical phase space meaningless on quantum level 
while configuration space retains its meaning. However, one can further demand that 
the coordinates are noncommuting operators which implies some uncertainty relations 
making also the notion of the point in space-time no longer sensible. The simplest way 
to do this is to impose the commutation rules 

where Q^" is a constant c-number tensor. Recently, there has been much activity 
concerning field theories on such noncommutative space-time [Tj, They appear 
to have some attractive properties. On the other hand their quantization seems to be 
more subtle problem than in the standard case. In fact, the noncommutative space-time 
can be replaced by its commutative counterpart provided one simultaneously replaces 
ordinary product of field variables by "star product" defined by 

^$2(3:) = e^®''"s^e^$i(a;)$2(2/) \x=y 

Therefore, once 7^ 0, the resulting Lagrangian contains time derivatives of arbitrary 
order; the theory is nonlocal in time. This makes the quantization procedure much more 
complicated. Indeed, within the standard framework, the first step to quantize a given 
classical theory is to put it in the Hamiltonian form. There exists the general algorithm 
which allows to construct the Hamiltonian formalism for higher-derivative jH], [l], (Hj 
and nonlocal |6], [7j, [8] theories. However, its main drawback is that the Hamiltonian is 
not bounded from below; the quantization can be formally carried out but the resulting 
theory has serious disadvantages like, for example, the nonexistence of stable ground 
(vacuum) state. This is the price one has to pay for the generality of Ostrogradski 
formalism. From this point of view it seems reasonable to pose the question whether, 
for a given specific system, there exist alternative canonical formalisms more adequate 
for quantization purposes. It can happen that, due to the peculiar properties of the 
system under consideration, there exists canonical formalism which, being quantized, 
produce quantum theory with more desirable properties than Ostrogradski approach. 
Our main motivation is to show, on the simplest example, that such a situation is 
possible; namely, that, in some cases, there exists a variety of Hamiltonians and the 
corresponding symplectic structures including those leading to the nice quantum theory 
(with stable ground state, etc.) 

The Ostrogradski instability is shared by all theories described by the Lagrangians 
containing time derivatives of at least second order. Moreover, the instability phe- 
nomenon seems to be not directly related to the nonlinear character of underlying 
dynamics. Therefore, the simplest model to be considered is the celebrated Pais- 
Uhlenbeck quartic oscillator |i9|, linear theory of fourth order. 
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We will study here alternative Hamiltonian formalisms for quartic oscillator. The 
starting point is the obvious observation that the general solution to Lagrangian equa- 
tion depends on four arbitrary constants. This implies that the corresponding Hamil- 
tonian system should have two degrees of freedom. By inspecting the explicit form 
of solutions wc find that there arc always at least two independent globally defined 
constants of motion which, in addition, are quadratic in dynamical variables. On the 
other hand, the Hamiltonian must be also a constant of motion. Therefore, we can write 
out the most general Ansatz for quadratic Hamiltonian. By demanding the canonical 
equations to be equivalent to the initial Lagrangian one we find the relevant Poisson 
structures. In principle, the family of candidates for Hamiltonian functions is much 
wider. First, one could take an arbitrary function of two above-mentioned constants of 
motion. Moreover,for some values of parameters the dynamics is superintcgrable, i.e. 
admits third independent globally defined constant of motion; then the most general 
Hamiltonian is a function of three integrals of motion. However, more complicated 
Hamiltonians result in more complicated or even singular symplcctic structures. This 
implies that the relation between basic dynamical variable, its time derivatives and 
Darboux coordinates is a complicated nonlinear one and it is not clear whether it can 
be promoted to quantum theory. 

Let us conclude this section with the following remark. Our construction is neither 
a pure application nor an extension of Ostrogradski algorithm. First, it cannot be 
applicable for all quartic systems. This can be seen by noting that we need here 
second ( independent of the Ostrogradski Hamiltonian) globally defined integral of 
motion. This implies that Ostrogradski dynamics is integrable which doesn't seem to 
be automatically true, in spite of the fact that the canonical equations for Ostrogradski 
Hamiltonian have a very specific form. Second, the family of Hamiltonians constructed 
here includes in some cases the positive-definite ones, the property not shared by Os- 
trogradski Hamiltonian. 

The details of our construction are presented in Sec. II while Sec. HI is devoted to con- 
cluding remarks. Appendix contains some additional remarks concerning the problem 
of embedding the fourth-order system into Lagrangian system of two degrees of free- 
dom. 



II Hamiltonian structures 



Our starting point is the following Lagrangian 




m .2 muj'^ 2 '^-^ •■' 






(2) 
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or, equivalent ly 



here 

2 1 ± VI - 4Aa;2 


The form of solution to eq. (jHI) depends on 2- There are the following possibili- 
ties: 

(i) < A < then 2 > ^1 7^ ^2! 

(ii) A = 0; the harmonic oscillator case 

(iii) A < 0; then ujI < 0, > 0, 0;^ + a;| 7^ 

(iv) A = then ul = uj2 = loj"^, i.e. we are dealing with degeneracy 

(v) A > both ujI ^ are complex 

2 1 ± zV4A^2 _ 1 
2 = > u;i = = Wo (5) 

We shall consider these cases separately. 



(1) The oscillatory regime ((i)) 

The general solution reads here 

q{f) = Aicos{uJit + cti) + A2Cos{uj2t + 02) (6) 



It depends on four arbitrary constants Ai^2, «i, 2 which can be found knowing 
q, q, q and Q at any given time. Consequently, there are at most four independent 
locally defined integrals of motion; however, at least one of them must depend explicitly 
on time. Two integrals can be readily found by computing 2 from eq. (0) and its 
first three time derivatives. In this way one obtains the global integrals (normalised 
for further convenience) 

Jl = -T^rr — ((^' +^1'^)' + ^'(^ + ^!<iy) 
= — i7 {(^ + ^liq + ^ky) (7) 

For generic values of parameters no additional independent globally defined integral 
(which does not depend explicitly on time) exists; our system is integrable but not 
superintegrable. However, for A, uP' such that ^ is rational, ^ = j, it becomes 
sup erintegr able. The additional integral can be constructed as follows [lOj. One writes 
sin{lai — ka2) = sin{l{u>it + ai) — k{u2t + a2)); the latter is expressible polynomially 
in sin{uji^ 2^ + «!, 2), cos{uJi^ 2^ + «!, 2) which, in turn, can be computed from eq. ^ 



4 



and its first three time derivatives. In what follows we are interested in generic values 
of A. Therefore, we consider Ji^ 2 to be the only relevant integrals. 

As usual, the integrals of motion are related to some symmetries. Using Noether 
theorem suitably generalised to higher- derivative theories one finds the symmetries 
responsible for the existence of Ji_ 2- They read 

q^q + e(q±{ujl-ul)q) (8) 

We can now construct the Hamiltonian formalism. There exists standard procedure 
called Ostrogradski formalism |lj which works for any higher-derivative theory. 
However, for a particular dynamics there can exist a variety of suitable Hamiltonian 
structures. 

For the reason explained in Sec. I we restrict ourselves to quadratic Hamiltonians. 
Keeping in mind that the Hamiltonian itself is a constant of motion and, moreover, its 
rescaling is equivalent to the time rescaling, one can write the following Ansatz 

H{I3) = Jicos/3 + J2sin(5, -n < f3 < n, (9) 

Using 

gW = {g("-i), H}, n = l, 2, 3 (10) 
one finds the following one-parameter family of Poisson structures 

[Q: Q} = 71 n + 



cos (3 sin (3 
{q, q} = 

^■} = -7(^ + ^) (11) 
cosp sm p 

2 3 
U2 uji 



cosp sm p 

{q, = 

cosp sm p 



with 



Let us note the following: 

(a) The Poisson structure exists for all j3 except 13 = —n, — |, 0, | (this can be easily 
understood from eq. both Ai A2 are needed to characterize fully the motion); 
consequently, there exists four disjoint sectors for f3 : (— tt, — |), (— f , 0), (0, |) and 
(^, tt). However, the transformation H —H, qi ^ pi leaves Hamiltonian equations 
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invariant. Therefore, it is sufficient to consider the sectors (— f , 0) and (0, ^). 
(b) One easily checks that 



det 



1 



(n)|3 



m, n=0 



cos(3 sinP 



(13) 



which is nonvanishing. We conclude that the Poisson structures fjll|) are sympletic. 
(c) For any admissible /3 



{ij + ujlq, q + uj2q} 

{'q +ujlq, q W2q} 

{q + uj^q, Q +ujlq} 

{q + ^Iq. q +^lq} 










(14) 



(d) The structures corresponding to different are different (i.e. not canonically 
equivalent). Indeed, q,q,q,and q are well-defined functions of canonical variables. 
Therefore, the canonical transformations cannot change the numerical values of the 
Poisson brackets. On the other hand, due to tuf ^ uj\, sin(3 and cos/3 are uniquely 
fixed once the RHS of eq. (Illj) are known. 

The canonical variables are found by passing to Darboux coordinates. There is a 
freedom in defining such a transformation - one can always perform an additional 
symplectic ( in our case - also linear ) transformation. We shall impose a further 
constraint Pi ~ q'j, i = 1,2. Using (c) one finds the following canonical variables: 
- for the (0, |) sector: 



qi 



6^ cos (3{q + ulq) 
Pi = m5\J cos P(q +ujfq) 
q2 

q - 



5^Jsml3{q + ujlq) 



m5\Jsm(3(q +ujlq) 
1 , qi 



(15) 



= ( 
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Hid) = (— H -qi) + (— H -qi) 



S = 



V2X 



, ,2 
UJ2 



-for the (— 1,0) sector: 



qi = 5\Jcosl3{q + ulq) 
Pi = m5\J cosl3(q +ciJ^g) 
52 = 5\J-sinl3{q + uolq) 
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P2 — —m5^—sinl3(<i -\-ujlq) (16) 

1 f qi ?2 X 



The formulae (15), (16) have a nice interpretation. The q variable is one of the coordi- 
nates of twodimensional quadratic system for which qi and q2 are normal coordinates 
[11]. However, for P G (— 1,0) the energy of one of the normal oscillations enters with 
negative sign. Actually, in each sector all systems look the same except the formula 
for q in terms of normal coordinates gi_2- 

Let us consider in some detail the sector (— |, 0). Passing to the Lagrangian 

m = {^qI - -^q!) - {^qI - -^qI) (17) 

is a regular procedure. On the other hand, under the canonical transformation 



^2 = , , „ ^ iPl + P2) (18) 




= ^\ 2 (^2^1 - ^1^2) 



the Hamiltonian attains the Ostrogradski form 

/f(/3)=P:,.-|^-y9l + — 5? (19) 

Note that eq. (19) defines now a singular Hamiltonian in the sense that momenta are 
not expressible in terms of velocities and coordinates; in fact, ?i = imply q2 = 

but also qi = q2, so pi cannot be expressed in terms of (ji, q^. Therefore, some care is 
needed when passing to the Lagrangian formalism which results in additional variable 
- the Lagrange multiplier enforcing gi = ^2; as a result qi obeys eq.(2). However, 
qi coincides with q only for [3 = — |; for other P's q and gi, are different linear 

combinations of normal coordinates gi,2- 

Using eqs. (16) one can express the Lagrangian (17) in terms of q variable (up to a 
total derivative) 

mS'^ 2 
L{f3) = —-{cosp + sin(3) q + 
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-m5'^{uJl{cos(3 + -sinP) + u;l{sinP + -cosP))(f + 
— ^((^1 + 2ujI)ujIcos(3 + (^2 + 2ujl)ujlsin(3)q^ + 

22/2 n , 2 ■ n\ 2 



UiUJ2{uJiCosf3 + uj2sin/3)q (20) 
which leads to the following equation of motion 

^2 ^2 ^2 

{{cosP + sin(3)— + {ulcosp + ulsinP)){— + '^?)(^ + ^Dq = (21) 

We see from eq.(j2H) that there appears a new mode a;^ = '^^'Zlp'^'^f^fi^'^ unless /3 = — |. 
This is not surprising. First, let us stress that the theory defined by eqs. (fT^ and (fTTj) 
solves the problem of finding the Hamiltonian system containing eq. (0) ( or, equiva- 
lently, eq. Q) as one of dynamical equations. Indeed, eq.© is the direct cosequence 
of the definition of q in terms of gi_2 and the basic dynamical equations the latter obey. 
Moreover, due to ui ^^1^2, in order to determine time-dependence of q one has to know 
both qi and q2 which implies one has to impose four initial conditions on q; therefore, 
the theory given by eqs. (fT^ and (fT7|) describes the general solution to eq.© for ar- 
bitrary value of f3. In fact, the present formulation does not differ very much in spirit 
from Ostrogradski formalism. In the latter one of the canonical equations implies that 
the substitution qi ^ q, q2 ^ Q is cosistent while in the former Hamiltonian equa- 
tions imply the consistency of the rule gi,2 ~ ^ + ^^^129- this reason, making the 
substitution ~ g + ufq in flT7|) one gets consistent equation in spite of the fact that 
the number of independent variables is reduced. On the other hand this is not a point 
transformation which in general results in new modes (see Appendix). 
Finally, let us compare our findings with those of Ref . [9] . Again, it is a matter of simple 
computation to verify that the formalism developed in Sec. IIA of [9] corresponds to 

(2) The case A < ((iii)) 

Let us pass to the case (iii). Putting Lof = — \ tui l"^ one gets 

(^-|^.P)(^+-.^k^O (22) 

with the general solution 

q = Ael^^l* + A'e"l^^l* + Bcos{uj2t + p) (23) 
As in the previous case one easily finds two integrals. 

h = — -M - I ^1 1' (if + ^2(9- I ^1 1' q?) 

-{(^+^lqf-\^i? {q + ^lqf) (24) 



m 



-U)o 



V2{\ |4 -ul) 
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However, there exists also the third globally defined integral. The reason for that 
is that there is now only one angle variable which has to be cyclic. The additional 
integral can be found by computing /n(e''^^l*) =\ uji \ t + InA and cos{uj2t + (3). Then 
C = arccos{cos{uj2t + (3)) — ■^^ln{Ae\'^'^\*') is time - independent and cosC can be 

computed from g, g, g and The resulting expression is rather complicated and will 

not be considered here. 

We proceed along the same lines as in the first case. Define the Hamiltonian 

H{(3) = hcos(3 + hsinP, -n < /3 < n (25) 
The family of admissible Poisson structures reads 



{q,q} = 7( ^ + —-5] 

cosp smp 



{9, = 7(- 



sin(3 cos [3 

m}--ii}-^-^) (26) 

smp cosj3 

{q, '1} = 

I 14 4 

7 = 



sinfi cosP 
1 



These structures can be obtained from eq.(ll) making the replacement o;^ — > — | c^i p. 
Again we have four sectors f3 and it is sufficient to consider two of them only: 
- for /? e (0, 1) one gets 



qi = 5^cosl3{q- I Ui ^ q) 
Pi = m5\J cosP(q — | p q) 
q2 = sin/3 {q + Jlq) 

P2 = mS^sinf3(q +0;^^^ ^27) 
^ ^-v^A(|a;i 

^^^^ = ^2;;^ + ^'^^ ^^2^- ^^^^^ 



5 = 

for f3 e (— 1,0) one gets 



-V2X 



qi = 6Jcos(3{q- \uJi\ q) 



Pi — m5\j cosP(q — I a;i q) 
q2 = 5^—sin(3{q + ulq) 

P2 — —m5^—sinl3((l +ujlq) (2^ 



^^^^ - ^2^^! + ^^^^ " ^2^^! 2^^^^ 

Again the conclusion is that the q variable is a linear combination of normal coor- 
dinates for some quadratic system. The only difference as compared with the previous 
case is that the forces are in part repelling. 

(3) The degenerate case ((iv)) 

Consider the double frequency case (iv): 

+ 2^^)^g = (29) 

Then 

q{t) = Aicos{V2ujt + Q!i) + A2tcos{\/2ujt + 0:2) (30) 
The relevant integrals of motion are (again suitably normalised) 

/, = ^((g +2u;2g)2 + 2uj\q + 2u^qf) 

/2 = ^(2(g+2a;'g)g-g' + 4a;V) (31) 

Again, our system admits third integral which is globally defined but complicated and 
won't be considered. 
We put 

H{(5) = hcosP + hsinP - % < p < n (32) 
, , cos(3 



and find 



2msm^/3 
{q,q} = 

{2cos(5 + sin(3)uj'^ 



{9,9} 



2msi'n?(5 



{2cosP + sinp)u;^ 

^^'^> = 2msin^P 

{q. = 

2{cosj3 + sinP)u!^ 
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There are now two sectors, (— vr, 0) and (0, n) and its sufficient to consider only one, 
say /3 e (0, tt). Note that for /3 = |, H = J2; this is possible because J2 depends on 
both Ai and A2. 

Let us define new variables 



\/sin(3 . 2 N 
qi = ^[q + 2u q) 

q^ = ^ J q + 2(1 + tgf^yq) (34) 



mcosd .... ^. ^. 9 ., 
m\fsvn3 ^ 



Then the Hamiltonian takes the form 



H^^ + mu;'{2qiq2-ql) (35) 
m 

while q is the linear combination of basic variables 

^ ^ cosP / qi ^ 

\ sin(3 cosjS ) 

The Hamiltonian does not depend explicitly on (5 and the only /^-dependence comes 
from the expression for q in terms of basic variables gi, 52- 

Complex frequencies ((v)) 

Finally, let us consider the complex frequencies case. Formally, one can use the 
results of (i) and define the integrals 

-^1 = r -A & +^0^)' + + '^hf) (37) 



V2(a;4 



■^2 = -^kT^-^M +^Uf + ^o(^ + ^0^)') 

They are no longer real but rather obey 

J2 = -Ji (38) 
The one-parameter Ansatz for the real Hamiltonian reads 

H{I3) = i{e"^J^ + e-^^Js) (39) 
The relevant Poisson structure is given by 

{g, q} = -i^{e^^ + e-^^) 
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{q, q} = 

{q, q} = t^{ule'P + c^o'e"*^) (40) 
{q. q} = -ilie'^uj^ + e-'^ul) 

{q. q} = 

{q, q} = -i^^ute^ + ute-'^) 
1 

Now, all values of — tt < /? < tt are admissible. Again, we could consider only half 
of this domain, say < /3 < tt, but there is no point to do this as we are dealing with 
one sector only. Define 

qi = ^{q + ^h) 

q2 = ^(^' + ^0) 

Pi = +^oq) (41) 

m .... _2., 
P2 = —{q +uJoq) 



e 



The Hamiltonian takes the form 



while the expression for q reads 

£qi-m2 ..ON 

? = -J^—f (43) 

The canonical variables are not real. In fact, qi = q2, Pi = P2- The real canonical 
variables are obtained by taking the real and imaginary parts 

91 = ■^{Qi + iQ2) 

q, = l=(Q,-tQ,) (44) 
Pl = ^iPl-^P2) 

P2 = ^{Pl+lP2) 

Then 

- (S + - (5 + !!M±^e5) + '^(^l-M. (45) 

Further change of variables transforming H{I3) into the sum of dilatation and rota- 
tion is also possible [Q, (cf. Sec. III). 
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Ill Concluding remarks 



Let us summarize our results. We have found essentially one-parameter families of 
inequivalent quadratic Hamiltonian structures in all cases (i), (iii) ^ (v). In the first 
two cases these families consist of four disjoint sectors while there are only two sectors 
in the (iv) case and one - in the (v) case. In each sector the Hamiltonian can be put 
into the parameter-independent form; the structures belonging to any sector differ in 
the way the g-variable is expressed in terms of basic variables. 

Due to the symmetry H ^ —H, qi ^ pi, one can reduce by two the number of 
sectors we have to consider. Therefore, in the case (i) one has basically two sectors. 
The Hamiltonian is, respectively, the sum or difference of two independent harmonic 
oscillators. Our q variable is a linear combination of two basic coordinates qi, q2- 
Taking into account the possibility of rescaling the Hamiltonian and performing simple 
canonical transformation qi ±qi, pi —>■ ±pi one concludes from eq. (fTH|) . (jlfij) 
that q can be arbitrary linear combination of qi, q2 except that both coefficients are 
nonvanishing. 

Similar results hold for the case (iii). The only difference is that now one oscillator 
describes the repelling linear force. In the degenerate case (iv) there is essentially 
one sector (if one again takes into account the symmetry qi ^ Pi, H ^ The 
Hamiltonian takes less familiar form ()35|) while q is given by ()36|). 

In the complex case (v) there is one sector even without taking into account the 
above-mentioned symmetry. The Hamiltonian is now the difference of two har- 
monic oscillators coupled by the interaction term proportional to the product of coor- 
dinate variables. Due to the fact that the kinetic energy is not positive definite passing 
to normal coordinates is now impossible. 

Obviously, the Hamiltonian structures considered in Ref. [H] are the particular 
elements of our families {(3 = for (i) and (iii), /5 = | for (iv) and [3 = — tt for (v)). 

Finally, note that, apart from the Hamiltonian, there is always an additional 
quadratic integral of motion. Therefore, we expect in all cases the separation of vari- 
ables is possible. This is obvious for the first two families. In the degenerate case one 
defines P 













= Ql P2 




muj 


Pi 


= p[ — rrujjq'2 


P2 


= 'mujq'2 



Then i/, eq. ()35p takes the form 

H = -ojiq[p'2 - q'2p[) - mu'iq[' + q'i) (47) 

which separates in polar coordinates. 

Finally, consider the case of two complex conjugated frequencies squared. Making 
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an Ansatz [H] 



1 



{{Qi -P2) - i{Pi - ^2)) 



Qi = 



1 



q2 = 



{{Qi - P2) + i{Pi - 0.2)) 



Pi = ^V^{{q2+pi) - Kp2 + qi)) 

P2 = ^v^((^2 + Pi) + i{P2 + qi)) 



(48) 



one obtains 



H = -( ){qip2 - q2Pi) + -(^0 - ^o){qiPi + 52^2) 



(49) 



i.e. the Hamiltonian becomes a commuting sum of angular momentum and dilatation 
and separates in polar coordinates. 

The Hamiltonian formalism provides the first step toward quantization. The stan- 
dard approach based on Ostrogradski formalism and Dirac procedure provides 
a consistent quantum theory. However, its serious drawback is that the quantum 
Hamiltonian is unbounded from below. One is not surprised that the Hamiltonian 
is unbounded from below if the classical motion is unbounded (the cases (iii) -7- (v) 
above). On the contrary, in the case (i) the motion is bounded while the Ostrogradski 
Hamiltonian is again unbounded. We have shown that in this case there exists the 
whole family of Hamiltonians which, after quantization, yield stable ground state. 
We have to stress that in all cases under consideration the quantization procedure is 
quite simple because the Hamiltonians are built with the help of operators well-known 
from ordinary quantum mechanics: oscillator Hamiltonian, angular momentum, dilata- 
tion operator etc. 

IV Appendix. 

Let us discuss in more detail the problem of embedding the fourth-order dynamical 
system into Lagrangian system with two degrees of freedom. First let us note the 
following. Assume we have the first-order Lagrangian 



L = L{q, q) 



(50) 



Let us make the following substitution 



q = q{x, X, x) 

L{x^ X, X, x) = L{^q{x^ x, x), q{x^ x, x, x)) 



(51) 



Then one easily derives the following identity 
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/ dq d dq dq\ f dL d ,dL \ 

We see that, in general, the new equation of motion contains additional solutions except 
the case ll = 0, If = 0; in the latter case (|51|) describes point transformation leading 
to the equivalent dynamics. 

Consider now the system of two decoupled degrees of freedom, 

L = Li(gi,gi) + L2(g2,^2) (53) 

Assume that 

qi = qiiq,q,q), i = 1,2 (54) 

be the substitution, in terms of one variable, consistent with the equations of motion. 
By the latter we mean that substituting into both equations 



dLj d ( dLi 



0, i=l,2 (55) 



a,F(g,(/,^-,g,g(^^)) (56) 



dqi dt \ dqi 
produce the same equation for q: 

dU_±(dLi\ ^ 
dqi dt \ dqi ) 

with some constans ai^2- Then implies for the Lagrangian 

L{(1, <i, q, q) = Li{qi, qi) + L2(g2, ^2) (57) 

the following identity 

dL_d^ fdL\ (P_ fdL\ _ d^ fdL\ _ 

( dqi d dqi d^ dqi\ . .. ,jy. 

^ S ""'dm^ "'df^^ J '''''' ' ^ ^''^ 

If it happens that a'^s are such that the second and third term on RHS of (j^Hj) vanish, 
L gives no additional solutions. 
In our case 

L = a, - ^) + «^ - ^) (59) 

and the consistent substitution reads 

gi ~ g + uj^q 
g2 ~ g + t^i? 

Then q can be expressed in terms of qi and q2 and one obtains the consistent embedding 
of q into two dimensional system of first order. Moreover, for ai = —02 L gives no 
additional mode. However, the procedure is consistent for any ai^2 provided ai ■ 02 7^ 
(we must have two degrees of freedom in order to be able to express q algebraically in 
terms of basic dynamical variables). 
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